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Introduction 

Many Earth-orbiting satellites contain instru- 
ments that are designed to provide measurements 
from which the vertical distribution of trace con- 
stituents in the Earth’s atmosphere may be inferred. 
Many of these (e.g., SAM I and II, SAGE I and II) ac- 
complish this by using the Sun as a source of radiant 
energy and measure the attenuation of this radiation 
as the Sun rises and sets with respect to the space- 
craft. The line of sight between the spacecraft and 
the Sun is referred to as the tangent ray. There is 
a point along the tangent ray where the height of 
the ray above the Earth’s surface is a minimum — 
the tangent point. This point, projected vertically 
to the Earth’s surface, locates the subtangent point. 
Most of the attenuation of the solar radiation occurs 
in the immediate vicinity of the tangent point, and 
hence the computed number density or mass density 
of the absorbing material is generally associated ge- 
ographically with the subtangent point and at the 
altitude of the tangent point. 

For most satellite geometries, during one com- 
plete orbital revolution there is one sunrise and one 
sunset event. (There are no sunrise and sunset events 
for some geometries in which the orbital plane of 
the satellite is nearly normal to the Earth-Sun line. 
These geometries are not considered further in the 
present paper.) A typical near-Earth satellite com- 
pletes 14 to 15 orbits per day, and so there are nomi- 
nally some 10000 sunrise and sunset events per year. 
The mission designer is concerned with the global 
distribution of the subtangent points locating these 
events, and for some types of missions wants to cover 
as much of the Earth’s surface as possible. Unfortu- 
nately, other mission constraints may be, and usually 
are, imposed on the mission, thus frequently present- 
ing the mission designer with a nonoptimal choice of 
available design parameters. 

The present paper describes a feasibility study 
that considers the use of stars as radiation energy 
sources to provide more flexibility to the design of 
occultation missions for atmospheric probing. 

The problem to be addressed in the present study 
is to investigate the answer to one basic question: 
Do enough stellar photons reach the detector of the 
instrument to produce a measurable signal propor- 
tional to the incoming photon flux? In order to an- 
swer this question, three basic parameters of the sys- 
tem must be examined: 

1. What is the number of photons/ (cm 2 -sec) 
reaching the collector optics of the radiometer sys- 
tem? This number is independent of the radiometer 
design and is a function only of the surface temper- 
ature of the star, its size, and its distance from the 


Earth. Assuming 100-percent efficiency in the optics 
and detector systems, this is the maximum number 
of photons available to produce a measurable signal. 
In practice, of course, only a very small fraction of 
these actually produce signal output. 

2. The efficiency of the optical system is not 100 
percent. So, of all the photons that enter the collector 
area of the optical system, how many actually reach 
the detector? This number is a function of the length 
of the optical path between the collector entrance and 
the detector, and of the number of times the optical 
path changes direction or is otherwise interfered with. 
The largest contributor to the loss of stellar photons 
is the type and number of filters placed in the optical 
path to isolate a particular wavelength region. The 
average collector efficiency is usually of the order of 
a few percent. 

3. Of all the photons that finally reach the detec- 
tor, how many can actually be detected and corre- 
lated with the number entering the collector area and 
hence related to the incoming radiative flux? This 
number is primarily a function of the type of detec- 
tor and of the quantum efficiency of the detector in 
the wavelength region under study. 

Item 1, being instrument independent, is the only 
part of this three-phased question that can easily 
be addressed and is the main thrust of the present 
paper. 

Items 2 and 3 are highly dependent on the phys- 
ical characteristics of a specific collector and detec- 
tor system and, consequently, can only be addressed 
meaningfully if reasonable numbers are known or as- 
sumed for its performance characteristics. 

Analysis 


If we assume that a given star radiates as a black- 
body and that certain astronomical parameters of the 
star are known, it is possible to compute the photon 
flux of the star outside the Earth’s atmosphere. 

The spectral flux emitted by the star at wave- 
length A, cm, and surface temperature T*, K, is (see, 
for example, Liou 1980) 


F\ (T.) 


27r/ic 2 1 erg 

A 5 exp (^fr-) - 1 cm 2 -sec-cm 


( 1 ) 


where 

h Planck’s constant, 6.625 X 10" 27 erg-sec 

c speed of light, 2.9979 X 10 10 cm/sec 

k Boltzmann’s constant, 1.38 x 10 -16 
erg/K 


The energy per photon of this radiation is he/ A, and 
thus if we divide equation (1) by this quantity, we 
get an expression for the photon flux at the surface 
of the star 


w.| = ^ 


photons 


A4 ex P iuk) ~ 1 cm2 - sec " cm 


( 2 ) 


This expression gives the total number of 
photons/sec emitted per unit area of the star, per 
unit wavelength, at the surface of the star. If 
R* is the radius of the star, the total number of 
photons/sec emitted from the star in all 
directions is 


N Tx (7Y) = 47 TRl N X (T , *) 


photons 

sec-cm 


(3) 


If we assume no attenuation of photons in the 
interstellar medium, the total number of photons/sec 
given by equation (3) passing through the sphere of 
radius R* must be the same passing through any 
sphere centered at the star. In particular, if d* is 
the distance of the star from the Earth, then the total 
number of photons/sec that were emitted by the star, 
and that cross a unit area at the Earth distance is 


N , * (T*) = 


Nt x (7U 
4n dl 


- (zf N ^ 


photons 
cm 2 -sec-cm 


(4) 


be measured directly in any of several wavelength 
regions — the magnitudes of importance to us are the 
visual and photographic magnitudes. 

The human eye and photographic film have dif- 
ferent spectral response functions that peak at dif- 
ferent parts of the visible spectrum. The eye is most 
responsive to radiation around 5290 A, whereas the 
spectral response function of photographic film peaks 
somewhat more toward the blue at about 4250 A. If 
we assume that the spectral response functions are 
flat, then the temperature of the star can be inferred 
from the visual and photographic magnitudes as fol- 
lows (much of the remaining development was taken 
in bits and pieces from essentially two references: 
Payne-Gaposchkin (1954) and Jeans (1961)— mostly 
from Jeans): From equation (1) the total energy ra- 
diated by the star in the wavelength interval A 2 — Ai 
is 


Ex — 47 tR% 


j* x 2 2t rhe 2 
Ai A 5 



dA ( 6 ) 

sec 


Near the middle of the visible radiation band, say 
5290 A, and for reasonable stellar surface tempera- 


tures of 3000 K to 16000 K 


exp 



» 1 


Let 0(A) be the spectral response function of some 
instrument (or of the human eye) with which we want 
to measure the energy from the star. Typically, 0(A) 
is defined for some definite wavelength region AA. 
It generally has a single maximum somewhere near 
the center of the band and drops off to zero as the 
edges of the band are approached and is defined to be 
zero outside the band. With this function, the total 
number of photons/sec to which the instrument is 
potentially responsive is given by 

»,<?.) = (r)T* A > ^ 

\ a* / J 0 cm z -sec 

(5) 

The major problem in evaluating equation (5) is, 
of course, that we seldom know the radius of the star 
R * and may not know its distance from the Earth 
d*. However, the ratio R*/d* can be estimated if we 
assume blackbody radiation and if we know two other 
parameters of the star: its visual magnitude and its 
surface temperature. The magnitude of a star can 


If we invoke this assumption, and use the mean 
value theorem of integral calculus, we can write 
equation (6) as 


Ex ~ 47ri?* 



(A 2 - Ai) 


erg 

sec 


(7) 


where A 0 is some wavelength, Ai < A 0 < A 2 . 

Stellar magnitudes are defined such that any- 
where along the brightness scale, a difference of 5 
magnitudes corresponds to a ratio of 100 in bright- 
ness, or 


mi — m 2 = 2.5 log 10 



( 8 ) 


where B here refers to brightness. For a given star, 
then 


m = -2.5 log 10 B + Constant 

— —2.5 log 10 Ex + Constant (9) 
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We can write from equation (7) 

l°g 10 e t = 2 log 10 R* 

- iog 10 [ ex P ( x~£rv )] + Constant 
0 625 

= 2 log 10 R* ~ ' 7 + Constant (10) 

Ao-/* 

and putting equation (10) into equation (9) 

1.5625 

m = -5 log io R* — — ; — h Constant (11) 

Aq 1 * 


For the Sun, m s = "26.72, d s = 1.496 x 10 8 km, 
and hence M s = 4.85. 

Since the magnitude definition, equation (8), ap- 
plies to absolute as well as visual magnitudes, we can 
interpret equation (11) as referring to absolute mag- 
nitude with, of course, a different constant. If we 
write equation (11) for a star and for the Sun and 
subtract to eliminate the constant, then 

1.5625 

M - M S = - 5 log 1Q R* + 7 rj, 

Ao 1 * 

r . „ 1.5625 

+ 5 log 10 R s ~ r = 

Aols 

Evaluate this at 5290 A and rearrange to obtain 


Now, visual radiation centers around 5290 A, and 
photographic radiation centers around 4250 A. If we 
put these values into equation (11) for A 0 , then the 
color index /, can be written 

7228 

1 = «*ph ot - m vis = — — + Constant (12) 

i * 


The constant is evaluated as follows: by international 
agreement, the color index is defined to be zero for 
type AO stars, whose temperature is known by other 
means to be about 11200 K. This gives a value of 
-0.64 for the constant, and equation (12) becomes 


/ 


7228 

~t7 


0.64 


(13) 


Since both m p h ot and m v \ s can be measured, the sur- 
face temperature of the star T* can be inferred from 
equation (13) -with all the proper caveats concern- 
ing the assumptions about blackbody radiation, cen- 
tering of the visual and photographic wavelengths, 
etc. 

The absolute visual magnitude M of a star is 
defined to be the magnitude of the star as seen from 
a distance of 10 parsecs (1 parsec, D , is 3.083 x 10 13 
km). Absolute magnitude is defined for both visual 
and photographic magnitudes. 

Let m\ — M, in equation (8), be the absolute 
visual magnitude. Then 


B\ = B 2 



and substituting this into equation (8) 


M — m = 2.5 log 10 
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(14) 


(15) 


logic, (fj) = ^ - 0.2M - 0.01 (16) 

We now use equation (15) to eliminate the absolute 
visual magnitude and write equation (16) as 

ios, ° (19 = W - o2m + iogi ° (I) ■ 101 

By again using the known values for the Sun and 
rearranging terms, we get 

logl ° {it) = - °- 2m - 8 ' 656 ( 17 ) 

Finally, it might be somewhat more useful to use the 
color index rather than the temperature, as this is 
tabulated in many almanacs for many stars. So, by 
use of equation (13) we can put equation (17) in the 
form 



0.817/ - 0.2m - 8.130 (18) 


This equation gives us a way of estimating the ratio 
R*/d* for a given star, given its visual magnitude 
and color index and the assumption that it radiates 
as a blackbody. The color index and apparent vi- 
sual magnitude for some 1480 stars are given in the 
Astronomical Almanac. 

No pretense is made, of course, that equation (18) 
will give the actual R*/d* ratios very accurately, 
because of the several simplifying assumptions that 
went into its derivation. It should, though, be at least 
radiatively consistent with the measured data and 
hence should perhaps be interpreted as an equivalent 
or effective value. 

Measurements of stellar diameters are rather 
sparse and imprecise. Measurements of stellar 
distance are much more accurate and consistent. 
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Table I gives these parameters for a few stars (Payne- 
Gaposchkin 1954), whereas table II uses these data to 
demonstrate the accuracy of equations (17) and (18). 
By comparing the last two columns of table II — i.e., 
the “measured” R*/d* with the calculated values — it 
is seen that equation (17) predicts this ratio surpris- 
ingly accurately, with only two glaring exceptions: 
/? Centauri and 61 Cygni A. Since the agreement is 
so good for the other stars, it can probably be consid- 
ered that either the published R*/R s for these stars is 
considerably uncertain, since the parallax measure- 
ments are generally of high quality, or that black- 
body radiation does not well characterize the radi- 
ation from these particular stars. Since there are 
11 instances of good agreement compared with only 
2 bad cases, equations (17) and (18) will be used to 
evaluate in the remainder of the present work. 

The paper by Ridgeway et al. (1980) was brought to 
the author’s attention by a reviewer. These writers 
discuss several modem techniques for measuring the 
angular diameters of stars (lunar occultation, speckle 
interferometry, etc.) and give the results of a num- 
ber of sets of measurements for 22 stars. Kitchin 
(1984) discusses the theory behind these techniques 
in great detail. Computations on these stars using 
equation (18) showed generally the same agreement 
as presented in table II for Payne-Gaposchkin’s data. 


Detection of Photons 

The number of stellar photons per unit time and 
area that reach the collector component of an optical 
instrument in a given wavelength band is given by 


N r (n) 


(t)T NMdx 


photons 

cm 2 -sec 


then the total number of photons that reach the 
detector of the instrument is 

Nd(T .) = a c £ c (^) 2 J* 2 4>{ X) N X (T,) dX 

= A c Ec <t>( Ao) AT e (T.) (20) 


where 4>(X 0 ) is some mean value of the spectral re- 
sponse function within the bandpass of the instru- 
ment. Since the energy of each photon is he / A erg, 
the total input power reaching the detector is 


Pd x (T*) — N d (T*) — x 10 7 watts (21a) 
1 A 0 

and the detector output power is 


P Do (T*) = P D . (T*)q watts (21b) 

where q is the quantum efficiency of the detector. 
With A 0 given in microns 

P D . (T*) - 1.986 x 10' 19 Nl) ( T *} watts (22a) 


P D o (T„) = 1.986 x 10' 19 qN ^ T *^ watts (22b) 

A 0 

The detector output current is 

i D (n) = Ni)(T*)qe amperes (23) 

where e = 1.6 X 10“ 19 coulombs, the charge of each 
electron. 


and with the same assumptions as made earlier, this 
can be written 


N C (T+ 



2 TTC 

1 

Id, ) 

A* 

ex p(^&:) 1 J 


x (A2 ~ Ai) 


photons 

cm 2 -sec 


(19) 


If we define 

A ( collector area, cm 2 

E r optical efficiency of collector 


Numercial Results 

In order to produce some meaningful numerical 
results, both the characteristics of a bright star and 
a realistic optical detector system must be assumed. 
The 15 visually brightest stars were selected as en- 
ergy sources. The parameters of these stars needed 
for the present study are given in table III. Of these 
stars, Capella, Arcturus, and Aldebaran are classi- 
fied as giants, whereas Betelgeuse and Antares are 
super giants. The rest are main sequence stars. The 
temperature range of 3000 K to 18000 K encom- 
passes practically all the visible stars and so these 
15 stars probably present a reasonable cross section 
of the stars available for an experiment of the type 
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described earlier. The visual magnitudes and color 
indices were taken from the 1985 edition of the As- 
tronomical Almanac. The effective surface temper- 
atures were computed from equation (13) and the 
ratio R*/d* from equation (18). 

The detector channels, bandwidths, and quan- 
tum efficiencies assumed in the present study are 
presented in table IV. It must be properly empha- 
sized that this selection is, at least for the purposes of 
this paper, completely arbitrary and for illustrative 
purposes only. 

Table V through IX present the numerical data 
computed from the data of tables III and IV and the 
analysis tools developed in the present paper. 

Table V presents computed values of N C (T*) and 
N d{T*) / A C E C from equations (19) and (20), respec- 
tively, since these two parameters are identical for 
unit response function. Recall that N c is the num- 
ber of photons per unit time and area that reach the 
collector. This number, is of course, independent of 
any instrument and is a function only of the prop- 
erties of the star and of its distance from the Earth. 
On the other hand, Np has units of photons/sec and 
represents the actual number of photons per second 
impinging on the detector. The quantum efficiency 
of the detector then determines how many of these 
photons are detectable and/or are turned into some 
other measurable quantity which is proportional (one 
hopes linearly) to N c . 

Tables VI through IX are perhaps of more util- 
ity. Table VI shows the detector input power per 
unit collector area in each of the channels and for 
each of the 15 selected stars and is computed from 
equation (22a). Table VII shows the detector output 
power per unit collector area. Table VIII presents 
the same data as table VII, but each entry has been 
divided by the assumed channel NEP (noise equiva- 
lent power) for the simulated instrument. Thus, each 
entry of table VIII is essentially the signal-to-noise 
ratio per unit collector area of the instrument for 
the selected channels. Table IX contains the detec- 
tor output current per unit collector area. It can be 
seen that even for channels 3 and 4, which have the 
highest quantum efficiencies, the detector response 
is marginal at best. Only Sirius and perhaps Cano- 
pus are bright enough to yield reasonable detector 
powers. 

Also, it must be realized that these numbers 
are for radiation detected outside the Earth’s at- 
mosphere. For an occultation experiment designed 


to measure the altitude distribution of trace con- 
stituents, the stellar radiation is further attenuated 
as the star rises or sets with respect to the spacecraft. 
When viewing the star along a ray whose closest dis- 
tance from the Earth’s surface is near the tropopause 
(10-12 km), the photon counts (table V) can be fur- 
ther reduced by more than two orders of magnitude, 
with a corresponding reduction in detector power and 
signal-to-noise ratio. 

Comments 

It must once again be emphasized that the num- 
bers shown in tables VI through IX were computed 
for a hypothetical instrument system. 

The author is not an instrument specialist; he 
is only mildly conversant with the fundamentals of 
collecting and detecting photons and has essentially 
no knowledge of the actual performance character- 
istics nor of the operating problems associated with 
real detectors. From conversations with such spe- 
cialists, the figures chosen to define the system char- 
acteristics used herein are reasonably close to what 
might be expected in the specialized literature, but 
neither the performance characteristics nor the con- 
clusions drawn from their use here should be consid- 
ered definitive. Consequently, the low detector pow- 
ers and currents displayed in the tables, which seem 
disappointingly low to this writer, may not at all dis- 
may an electronic specialist who is working at or near 
the state of the art in photon detector systems. The 
use, for example, of photomultipliers in the visible 
portion of the spectrum, and avalanche or cascade 
photodiodes in the near infrared, may indeed make 
the detection and amplification of such low power 
levels possible and useful. The author thus draws 
no feasibility conclusions from the tables and leaves 
their interpretation and utility to those more famil- 
iar with operational detection systems. Table V is 
presented without prejudice; it certainly defines the 
upper limit in what can and cannot be detected and 
presents initial design information for the design of a 
working instrument. 


NASA Langley Research Center 
Hampton, VA 23665-5225 
April 21, 1987 
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Table I. Selected Parameters for a Few Stars 


Star 

m 

M 

T*,K 

R* / Rs 

d*/D 

f3 Centauri 

0.9 

-3.8 

21000 

ii 

0.762 

Sirius A 

-1.6 

1.3 

11200 

1.8 

.376 

Procyon A 

.5 

3.0 

6 500 

1.9 

.291 

a Centauri A 

.3 

4.7 

6000 

1.0 

.756 

61 Cygni A 

5.0 

8.4 

3 800 

.7 

.296 

Kruger 60 A 

9.2 

11.2 

3 300 

.34 

.256 

Barnard’s Star 

9.7 

13.4 

3100 

.16 

.543 

Arcturus 

.2 

-.3 

4 100 

30 

.080 

Aldebaran 

1.1 

-.1 

3300 

60 

.057 

(3 Pegasi 

2.6 

-1.4 

2 900 

170 

.016 

Betelgeuse 

.9 

-2.9 

3100 

290 

.017 

Antares 

1.2 

-4.0 

3100 

480 

.0095 

Sirius B 

8.4 

11.2 

7 500 

.034 

.376 


Table II. Measured and Calculated i?*/d* Values for the Stars of Table I 


Star 

Z?*, km 

d 

*, km 

(/?*/ d*) meas 

(R*/d*) 

eq. 

(17) 

(3 Centauri 

7.66 

X 

10 6 

4.05 

X 

10 13 

1.89 

X 

10' 

-7 

2.79 

X 

io- 

-9 

Sirius A 

1.25 

X 

10 6 

8.20 

X 

10 13 

1.52 

X 

10* 

-8 

1.56 

X 

IO' 

-8 

Procyon A 

1.32 

X 

10 6 

1.06 

X 

10 14 

1.25 

X 

io- 

-8 

1.42 

X 

10' 

-8 

a Centauri A 

6.96 

X 

10 5 

4.08 

X 

10 13 

1.71 

X 

10' 

-8 

1.86 

X 

10- 

-8 

61 Cygni A 

4.87 

X 

10 5 

1.04 

X 

10 14 

4.68 

X 

io- 

-9 

7.92 

X 

10- 

-9 

Kruger 60 A 

2.37 

X 

10 5 

1.20 

X 

10 14 

1.97 

X 

IO' 

-9 

1.97 

X 

10- 

-9 

Barnard’s Star 

1.11 

X 

10 5 

5.68 

X 

10 13 

1.96 

X 

10' 

-9 

2.04 

X 

10- 

-9 

Arcturus 

2.09 

X 

10 7 

3.85 

X 

10 14 

5.42 

X 

io- 

-8 

5.56 

X 

10- 

-8 

Aldebaran 

4.14 

X 

10 7 

5.41 

X 

10 14 

7.73 

X 

io- 

-8 

8.20 

X 

10- 

-8 

/ 3 Pegasi 

1.18 

X 

10 8 

1.93 

X 

10 15 

6.12 

X 

10" 

-8 

7.26 

X 

10- 

'8 

Betelgeuse 

2.02 

X 

10 8 

1.81 

X 

10 15 

1.11 

X 

io- 

-7 

1.17 

X 

10- 

-1 

Antares 

3.34 

X 

10 8 

3.25 

X 

10 15 

1.03 

X 

io- 

-7 

1.02 

X 

10" 

~7 

Sirius B 

2.32 

X 

10 4 

8.20 

X 

10 13 

2.83 

X 

10“ 

-10 

2.83 

X 

10- 

- 10 



Table III. Visual Magnitude ra, Color Index /, Surface Temperature T*, Radius/Distance 
to Earth Ratio R±/d* for Sun and 15 Bright Stars Used in Text 


Star 

m 

I 

T*. K 

R*/d* 

Archernar 

0.46 

- 0.16 

15058.33 

0.44389474 X 10~ 8 

Aldebaran 

.85 

1.54 

3315.60 

.90819687 x 10~ 7 

Rigel 

.12 

-.03 

11849.18 

.66296408 x 10 -8 

Capella 

.08 

.80 

5019.44 

.32181034 x 10 -7 

Betelgeuse 

.50 

1.85 

2902.81 

.19118332 x 10“ 6 

Canopus 

-.72 

.15 

9149.37 

.13694620 x 10 -7 

Sirius 

- 1.46 

.00 

11293.75 

.14521116 x 10 -7 

Procyon 

.38 

.42 

6818.87 

.13713238 x 10“ 7 

Spica 

.97 

-.24 

18070.00 

.30193955 x 10~ 8 

Hadar 

.61 

-.23 

17629.27 

.36315330 x 10~ 8 

Arcturus 

-.04 

1.23 

3865.24 

.76367751 x 10 -7 

Rigel Kent 

.33 

1.00 

4407.32 

.41783037 x 10 -7 

Antares 

.96 

1.83 

2926.32 

.14897384 x 10 -6 

Vega 

.03 

.00 

11293.75 

.73113908 x 10 -8 

Altair 

.77 

.22 

8404.65 

.78657475 X 10 -8 

Sun 

- 26.72 

.55 

6073.95 

.46062765 x 10 -2 


Table IV. Midpoints of Detector Channels, Bandwidths, and Quantum Efficiencies 
Assumed (NEP = 1.4 x 10 -14 watts for Each Channel) 


Channel 

A 0 , //m 

AA, ^m 

q, electrons/photon 

1 

1.0197 

0.0196 

0.21 

2 

.9355 

.0200 

.35 

3 

.5999 

.0145 

.76 

4 

.5250 

.0148 

.65 

5 

.4524 

.0019 

.42 

6 

.4475 

.0032 

.41 

7 

.3846 

.0198 

.30 
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Table VIII. Detector Output Power per Unit Area Divided by Channel NEP (1.4 X 10 watts) of Detector Pjj q (T*)/A c E C NEP 
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